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ABSTRACT 

We present the first results of a survey of 14 low redsliift galaxy clusters 
using Suzaku. Although luminous {L^ > 1 x 10"^^ erg s~^ (0.1-2.4 keV)), these 
clusters have no prior pointed X-ray data. Together with 47 other systems they 
form a flux limited sample {fx > 1.0 x 10^^^ erg s^^cm"^ (0.1-2.4 keV)) with 
z < 0.1 in the northern celestial hemisphere. Using this total sample we evaluate 
the local L — T relationship and the local cluster temperature function. Suzaku 
temperature measurements appear to be in accord with those of other missions. 
General agreement is found with other published estimates of the low redshift 
cluster temperature function, however the sample used here exhibits slightly lower 
space densities at gas temperatures below 4-5 keV. We find a corresponding deficit 
in the number of clusters with temperatures between approximately 4.5 and 5.5 
keV. Although at low significance, a similar feature exists in previous low redshift 
cluster datasets. We suggest that low-redshift cluster samples, while crucial for 
calibrating precision cosmology measurements, must be used with caution due to 
their limited size and susceptibility to the effects of cosmic variance. 

Subject headings: cosmology: observations - X-rays: galaxies: clusters - inter- 
galactic medium 



1. Introduction 



The nature and evolution of galaxy clusters offers an indepen dent and complement ary 
probe of cosmology to those of the cosmic microwave background (IHinshaw et al. 20081 ) or 
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supernova (e.g. Permlmutter & Schmidt 2003). Theoretical work (e.g. Haiman, Mohr & 
Holder 2001) has shown that large surveys of clusters to reasonably high redshift {z ~ 1) 
can not only provide these much needed constraints, but can, through the use of information 
such as the spatial distribution of clusters, be made "self-calibrating" - if the evolution of 
cluster structure can be parameterized (Hu 2003, Majumdar & Mohr 2003), and if the scatter 
between cluster mass and observables (luminosity, temperature) is understood for clusters 
at the survey detection threshold. 

Future X-ray surveys or Sunyaev-Zel'dovich Effect (SZE) surveys, such as the South 
Pole Telescope (SPT), and the Atacama Cosmology Telescope (ACT), will be designed to 
detect ~ 10^ clusters to high redshift within a limited sky area and will not sample very 
many local clusters, or obtain very high physical resolution data. However, both of these 
characteristics are vital in order to establish the zero-redshift baseline from which we can 
then evaluate the global evolution of cluster scaling laws and structure, and calibrate possible 
systematic effects induced by potentially complex astrophysics. 

In this paper we describe the results of an effort to utilize existing X-ray cluster data, 
and new observations by the Suzaku (JAXA/NASA) mission to add to the existing body 
of work on the low redshift population of clusters (e.g. Henry 2000; Reiprich & Bohringer 
2002; Ikebe et al. 2002). Specifically, we present the results of a Suzaku survey of 14 nearby. 
X-ray luminous, clusters that have otherwise never been observed in a direct, pointed X-ray 
observation with earlier, spectrally sensitive instruments. Our motivation is to increase the 
robustness of the local cluster temperature function at lower masses, which are typical of the 
mass threshold of the majority of likely future X-ray and SZE cluster surveys. The range 
of masses probed by this data combined with pre-existing data (a total, complete, sample 
of 61 objects) is approximately 8 x 10^^ — 10^^ Mq {Hq = 70). We therefore bracket the 
expected low mass limit of future surveys, all of which are designed to probe to ~ 2 x 10^^ 
Mq (locally for the X-ray, and across essentially all redshifts for SPT/ ACT). For modern 
cosmological tests based on the form of '^'^^^^^^^ (the number of clusters per unit redshift 
above a mass M, c.f. Haiman, Mohr & Holder 2001), the only theoretical requirement is a 
precise understanding of the scatter in mass around the detection threshold (e.g. X-ray flux 
limit). 

The rest of this paper is organized as follows. In § |2l we compile the local, fiux-limited, 
complete sample. In § [3l we describe the Suzaku data, and the spectral analysis method. 
We also present spectral fitting result for the 14 clusters. 

In § m we fit the L — T relation in two different ways. In § [5l we describe our procedures 
for computing the XTF and present our results. We also make a comparison with previous 
results. In § [6l we briefly summarize our conclusions and the implications of this work. 
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Throughout this paper, we adopt a spatially flat universe dominated by a cosmological 
constant and cold dark matter (CDM), with the following set of cosmological parameters: 
flrn = 0.27, Qa = 0.73 and Hq = 70 km s'^ M-pcK 

2. Local Flux-limited Cluster Sample 

Our sample is derived from the ROSAT Brightest Cluster Sample (BCS) (Ebehng et al 
1998). The BCS is based on a pure X-ray flux selection criterion using the ROSAT all-sky 
survey data. In the northern celestial hemisphere the BCS is (in parallel with the REFLEX 
survey in the south, Boehringer et al. 2001) the most complete, flux limited sample of low 
redshift clusters to-date. 

In Figure 1 we plot the luminosity versus redshift for BCS clusters; the fiux hmit is 
4.4 X 10"^^ erg s"^ cm"^ (0.1-2.4 keV) and all objects have Galactic latitudes \b\ > 20°, and 
redshifts z ^ 0.3. For the local cluster sample proposed here we have chosen a subset of the 
BCS with a flux limit of 1 x 10^^^ erg s^^ cm~^ (0.1-2.4 keV) and a redshift limit of z = 0.1 
(shown in Figure 1). Making this fiux and redshift cut maximizes the number of objects 
that have prior X-ray pointed imaging spectroscopy data while minimizing the number of 
new datasets required to complete a statistically useful sample. For example, the flux limit 
required to ensure a 100% completeness in preexisting spectroscopic data is ~ 2.3 x 10~^^ 
and would result in a sample of only 20 clusters, with only 3 of these with < 10^^ erg s~^, 
compared to 18 below this luminosity in our sample. Our chosen redshift constraint ensures 
that our sample is in the local universe, minimizing any redshift dependent variation/scatter 
in the L-T relation and XTF which are to be calculated. We note that our fiux limit is slightly 
lower than the 2 x 10~^^ erg s~^ cm~^ limit of Ikcbc ct al (2002, hereafter 102), although 
their fluxes were recomputed from the ROSAT data and are therefore slightly different (see 
also Reiprich & Bohringer 2002). Our final sample contains a total of 61 objects (see below), 
14 of which we have observed in Suzaku pointed observations (§3). 



3. Suzaku Data and Analysis 

Suzaku is Japan's fifth X-ray astronomy mission, and is equipped with an X-ray Imaging 
Spectrometer (XIS) and hard X-ray detector (HXD). The 14 clusters identified as having no 
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existing pointed X-ray imaging spectroscopy data were observed using the Suzaku XIS over 
an approximately 12 month period from 2006-2007. Table 1 summarizes these observations. 

The Suzaku XIS consists of four co-aligned X-ray CCDs (Koyama et al. 2007). Three of 
these are front-illuminated (XIS0,2,3), the other is back-illuminated (XISl). The CCDs are 
mounted at the focus of four independent X-ray telescopes (XRT). Each CCD is 1024 x 1024 
pixels, and covers an image region of 18' x 18' on the sky. The PSF of the XRT has a half 
power diameter of 2'. Sensitivity is in the energy range of 0.2-12 keV. 

Exposure times were chosen to ensure that approximately 15,000 XIS counts in total 
(3 FI + 1 BI) were obtained in the cluster core region for each object. A ratio of flux in 
the core to the total flux (integrated to infinity for a /3 = 2/3 mean emission profile) is a 
constant 29%. The estimated core flux from the BCS was then further reduced by 50% to 
allow for vignetting of a diffuse source across the XIS f.o.v. 

Not all of the 14 clusters were observed by all 4 cameras, 5 clusters (Z8852, A2665, 
A2495, A2249 and A566) were not observed by XIS2. We show XISO (raw) images of the 14 
clusters in Fig. [2] together with brightness contours. The color and contour are both on a 
logarithmic scale, and the data are binned by a factor of 4, then smoothed by gaussian with 
a = 5 pixels. 

We start our spectral analysis from cleaned event files, using the standard cleaning 
process (selection of the event grade and bad CCD column, disposal of non-observational 
intervals, and removal of hot and fiickering pixels). We then accumulate the spectra within 
a radius of 7.2 arcmin of the nominal cluster center (the peak of surface brightness) using 
the "xselect" Ftool. The spectrum is rebinned by a factor of 8 so that there are sufficient 
events in every energy bin for robust fitting. Blank sky spectra are used as the background, 
and subtracted from the cluster spectrum. The net exposure time of blank sky spectra 
is 96.1 ks. The galactic foreground emission is ignored. This is justified by the fact that 
our clusters are all hot (T > 3 keV as we will see below), and the galactic latitudes of 
Suzaku clusters and blank sky are all over 20 degrees. Thus the amount of galactic emission 
over the relevant energy range (i.e. where the largest photon count is) is very small in 
comparison. We generate the Redistribution Matrix File (RMF) and Ancillary Response 
File (ARE) by using the tools "xisrmfgen" and "xissimarfgen" included in the HEASoft 
package (Ishisaki et al. 2007). The XIS image is used as the input surface brightness 
distribution for "xissimarfgen". The HEASoft package used is HEAsoft 6.3.0. Note that we 
use a fixed angular size for accumulating spectra, i.e. higher redshift clusters are sampled 
with larger physical radius. Recent observations indicate that cluster temperature profiles 
generally decline with radius outside the core (Vikhlinin et al. 2005, Leccardi & Molendi 
2008), this implies that in our analyses hotter regions are weighted preferentially for lower 
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redshift clusters. We therefore adopt a phenomenological cluster model that is a reasonable, 
albeit approximate, fit to currently observed temperature profiles in order to estimate this 
uncertainty (Reid & Spergel 2006, Shang et al. 2008 in preparation). Wc find that the 
mean emission weighted temperature within 7.2 arcmin of a cluster center and that outside 
the core (~ 0.05-0.1 virial radius) differ up to 4% for z=0.0341 (the lowest redshift in our 
sample) and z=0.0980 (highest redshift in our sample). The deviation of these temperatures 
from the mean temperature within R^qq is also smaller than 4%. While there may be an 
effect we therefore consider it to be smaller than other sources of experimental uncertainty 
- either systematic or random. 

Ideally the measured intracluster gas temperature would correspond to the virial tem- 
perature of cluster. However, clusters are not always isothermal. Merging, cooling, and 
phenomena such as AGN outflows can disrupt isothermality. Approximately half of all sys- 
tems exhibit evidence for being cooling core clusters (CCC), which have a measurably cooler 
gas component in their central regions (e.g. Peterson et al. 2003). The cooling core can 
contribute a significant fraction of the luminosity of cluster, so a simple emission-weighted 
mean temperature may be considerably lower than the true virial temperature. Excision of 
the central emission is non-trivial in the Suzaku data due to the PSF. A more appealing 
solution is to adopt a suitable spectral model. 

In our analyses, we test both single temperature (IT) and two temperature (2T) models 
as fits to the net spectrum. A IT model should provide a better fit to non-cooling core, 
isothermal clusters (NCCC), while a 2T model should offer an improved fit to a CCC. The 
2T model represents 2 independent gas components: a cool component from the core and a 
hotter component which should be more representative of the virial temperature. 

The MEKAL model is taken as our plasma code, as implemented in XSPEC version 
12.3. Ix. We simultaneously fit the spectra accumulated from the different CCD cameras. 
The energy range used is restricted to 0.5 keV to 7 keV to mitigate some of the contamination 
issues at lower energies and to minimize the background at higher energies relative to the 
cluster emission. We also tried ignoring the energy band around the Si K-edge (1.825-1.840 
keV) due to a potentially anomalous response in that band, however the difference in the 
resulting temperature is negligible (< 1%). The temperatures are tied together, whcras the 
normalizations are treated independently for different cameras (however, the ratios between 
the normalization of the hot component and that of the cool component in the 2T model 
are also tied). We try both fixing the metal abundance and setting it as a free parameter. 
When the abundance is fixed, the value is set to be 0.3 of solar. The galactic HI column 
density is fixed at the value obtained from the LAB survey (Kalberla et al. 2005). We also 
enforce an initial upper limit to the hot component :T/jot in our 2T model, since it is assumed 
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to represent the temperature of the majority of the ICM gas. The results of the spectral fits 
are summarized in Table [2] (metallicity fixed) and Table [3] (metallicity as a free parameter). 
Example spectra are also shown for clusters A1800 and Z8276, together with IT model fit 
(metallicity fixed) in Fig. [31 The fits are good overall, except for the low energy range (< 1 
keV) and some residual Fe features around 1 keV and 6.7 keV for Z8276. We suspect the 
mismatch below 1 keV is due to uncertainty in contamination in the optical blocking filter 
of Suzaku XIS (Sato et al. 2007). We estimate this uncertainty by fitting the spectrum 
using the energy range 1 keV-7 keV, and comparing the result with that from 0.5 keV - 7 
keV. For the majority of our clusters, the temperature difference is small (~ 2 — 3%), but 
for a few such as UGC03957 and A2495, the best fit temperature could drop by up to 7%. 
The remaining Fe feature in Z8276 appears to be because we fix the metallicity to 0.3 Zq, 
while the true value must be higher than this. This feature disappears when we thaw the 
metallicity. The resulting temperature is only slightly altered (see Table [2] and Table [3]). 
We also tried replacing the mekal model with the vmekal model, and allowing He and Fe 
abundance to vary, again the difference in temperature is very small (< 1%). 

To choose the most likely virial temperature of the cluster from our IT and 2T models, 
we perform an F-test that indicates which model is more applicable for the cluster by com- 
paring and the number of free parameters.. If including another temperature component 
(2T model) gives rise to a more than 99% significant improvement, the higher temperature 
(Thot) from the 2T model is chosen as the temperature for the cluster, which also implies 
it is a CCC Otherwise, the IT model temperature is picked as representing the virial tem- 
perature, and the cluster could be viewed as NCCC. However, there is one system: A2665, 
where although the 2T model passes the F-test, we still choose the result from IT model. In 
this case Thot from the 2T model for A2665 is more than 3 times the temperature obtained 
from the IT model fit, and the normalization of the hot component is much smaller than 
that of cool component. For a typical CCC, it should be larger or at least of the same order. 
So we do not believe that T^ot is actually representive of the virial temperature for A2665. 
The model chosen for the virial temperature is indicated in bold face in Table [2] and Table 
[3l The temperatures shown in Table [2] are used in our analyses below. 

Comparing the results in Table [2] and Table [31 we see that the temperature generally 
varies only a little after setting the metallicity as a free parameter instead of fixing it at 0.3 
solar. The one exception is UGC03957. This system is better fit as a NCCC after we thaw 
the metallicity parameter, it also has the highest metal abundance in our Suzaku survey 
(0.59 Zq). 

To check whether the Suzaku temperatures obtained using our procedure are consistent 
with published value from other missions, we also performed our analysis on 3 clusters 
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(A1060, A2218 and A1795) with well established temperatures in the literature from multiple 
missions (e.g. ASCA, Chandra, XMM). Our results appear to match previous results quite 
well in both NCCC and CCC cases. A detailed comparison is listed in Table HI Below are a 
few notes for each cluster. 

A1060: With z=0. 00114, it is one of the nearest Abell clusters. Various observations 
show that it has uniform temperature and abundance distribution, and is therefore thought 
to be in a dynamically relaxed state (Tamura et al. 1996, Furusho et al. 2001, Hayakawa et 
al. 2004, Hayakawa et al. 2006). Although recent observations indicate that the temperature 
profile has a steep drop in the outer regions (Sato et al. 2007), the variation within 7 arcmin 
is small. Obviously, a simple IT model is sufficient and perfectly suitable for such a case. 
All temperature measurements give quite close results: 3.26 it 0.06keV from ASCA(Furusho 
et al. 2001), 3.36^°;°^ keV from Chandra (Hayakawa et al. 2004), 3.27 ± 0.02 keV from 
XMM-Newton. Our result is also in perfect agreement with them. 

A2218: Residing at z=0.175, this cluster is further away than any object in our sam- 
ple. Various observations suggest it is at a later merger stage with a centrally peaked and 
azimuthally asymmetric temperature map (Govoni et al. 2004; Machacek et al. 2002; Marke- 
vitch et al. 2000). A IT model is preferred for this cluster since there is no cooling flow in 
the center, and our result is in close agreement with previous measurements. 

A1795: Is a strong cooling flow cluster at z=0.0634. Different measurements and dif- 
ferent treatments of the cooling flow give very different results, varying over a wide range 
of 5.3-7.8 keV. The low end of the measured temperature range (5.3keV) is from IT models 
that do not take into account of cooling flow. But even within the analyses that explic- 
itly considered the cooling flow (fltting with a 2T model, cooling flow model, or excising 
the coohng flow region), the range of results still appears to be large. Fortunately, recent 
results seem to stablilize around 6 keV. Among them, Allen et al. (2001) extracted the tem- 
perature from ASCA data using 5 different models. Their result is summarized as follows: 
5.40^Q09keV from a IT model (absorbing column density flxed), 5.33^Q}°keV from a IT 
model (absorbing column density as a free parameter), 5.87lo igkeV from a constant cooling 
flow model, 6.2llQ22keV from a 2T model, an d 5.6llQ'J5keV from a isothermal cooling flow 
model. An incomplete list of other measurements includes: 6.1keV from Chandra (Vikhlinin 
et al. 2005, IT model with central 70 kpc excluded), 5.95 ± O.lkeV from XMM-Newton 
(Arnaud et al. 2001, IT flt with cooling flow region excluded). Our result (5.73tQ;2okeV is 
in good agreement with these recent measurements. 

The ROSAT fluxes used to deflne our sample are computed from raw counts assuming 
a spectral model and gas temperature derived in an iterative fashion from an assumed L — T 
relation (Ebeling et al. 1998). We have therefore checked whether our spectroscopically 
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derived temperatures alter the cluster flux sufficiently in any case to cause systems to fall 
below our flux limit. We find that the fluxes vary only some 1-2 % even if the temperature 
changes by a factor of 2. Since this is small compared to the ROSAT count rate uncertainty 
(around 15%) it is unlikely to bias the sample either by excluding or including systems at 
the flux limit. The ROSAT flux is therefore shown in Table [5l but luminosities have been 
computed using our updated temperatures where appropriate. 

57 clusters in Table O have measured/published temperatures, which is 93% of our 
sample. 34 clusters were measured by ASCA, 14 were measured by Suzaku, 6 by Einstein 
or EXOSAT, 2 by Chandra, 1 by XMM-Newton, 4 used estimated temperature from our 
L — T relation which is to be calulated. When a cluster has temperature measurements 
by several instruments, our ffist choice is to adopt the ASCA measurement in order to 
minimize any instrumental systematics. More than half of the clusters in the sample have 
ASCA observations, in addition the ASCA SIS is very similar to the Suzaku XIS. Our second 
choice is to use Chandra or XMM-Newton for their accuracy, lastly we choose others. If the 
spectrum of a cluster was fitted with different models in the literature, we choose 2T model 
fits over cooling flow models to ensure consistency with our analysis. 

Our sample differs from that of the HIFLUGCS (Reiprich&Bohringer 2002) used by 
102 in the foUwing aspects: HIFLUGCS covers both northern and southern skies, ours only 
covers the northern sky; HIFLUGCS uses a flux limit of 2.0 x 10~^^erg s~^ cm~^, while ours 
is a factor of two lower; the two samples are constructed based on ROSAT flux estimated in 
slightly different ways. 

4. The L-T Relation 

The X-ray luminosity-temperature (L — T) relation for galaxy clusters has long been 
used as a probe of their physics and evolution. In this section, we fit our data above with 
the following relation, 

logLo^_2^1,gV(^"' ergs s-^) = a + /31ogT(keV) (1) 

The L—T relationship has an intrinsic scatter in addition to any measurement errors. A com- 
mon way of dealing with this kind of problem in curve fitting is to employ the BCES(X2|Xi) 
estimator (Akritas & Bershady 1996). However, the slope derived from the BCES estimator 
depends on the measurement error. In fact, the data with the largest errors/scatter are 
effectively given more weight in the fit than the data with smaller errors/scatter, contrary 
to naive expectations. In an attempt to mitigate this inherent bias we also performed a 
modified fit which takes intrinsic scatter into account explicitly. Our is defined as the 
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following: 



i logLi ' logTi mtr 



where a? , is the intrinsic scatter, and a? and a? are the measurement variances 
mtr ' logL, logTi 

in L and T. cr?^^^ is chosen so that the minimum reduced is 1- The results of both a 
standard BCES fit and that of our modified ^ire listed in Table [Hi 

As we can see, the two methods are roughly consistent for our data, but the fit h^is 
the advantage that small error measurements are given more weight in the fitting. This 
approach yields slightly lower power law indices for the L — T relation. 

In the analyses described below, we use the result from the fitting. The fits together 
with data are plotted in Fig. HI 

In the same way we also fit the relation between bolometric luminosity and temperature. 
These results are also presented in Table El the measured slope steepens to ~ 3. The 
conversion between band and bolometric luminosity is computed using the Mekal plasma 
code and measured temperature. Compared with high redshift measurements (e.g. Ota 
et al. 2006), the slopes match within the uncertainty, while our normalization is smaller, 
which is in general agreement with the positive redshift evolution of luminosity that has been 
suggested by previous works (Vikhlinin et al. 2002, Lumb et al. 2004, Kotov & Vikhlinin 
2005). 



5. X-ray Temperature Function 

5.1. Derivation of XTF 

As a first step towards using our cluster sample to provide a complementary zero-redshift 
calibration point for cosmological studies (alongside, for example, Reiprich & Bohringer 2002, 
Ikebe et al. 2002) we have made a preliminary assessment of the cluster X-ray temperature 
function (XTF). 

The number density of clusters hotter than a temperature T is given by, 

niT) = E (3) 

where Vmax{T) is the maximum comoving volume where a cluster of temperature T would be 
detected. Following many previous works, we assume that logL follows a normal distribution 
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with mean value given by the L — T relation above, and with variance cr^^j^^j.. 




(4) 



where the Vmax{T, L) is the maximum search volume for a cluster with temperature T and 
luminosity L. 



where the Q is the solid angle sky coverage of the survey (equal to 4.14 for the ROSAT 
BCS), and ^max = Min{zf, Zs). Zg is the redshift limit of our sample, which is 0.1. Zf is 
constrained by the flux limit, and is computed by solving the following equation. 



where L^^^{L,T) is the observed luminosity after applying a K-correction, and f\[^ is the 
flux limit of our sample, which is 1.0 x 10~^^erg cm~^ s~^. Since we assume a fiat universe 
model, the comoving distance takes a simple form. 



The calculated VmaxiTi) is listed in Table [51 We plot the cumulative XTF in Fig. 
[5] together with the results from previous works. We find general agreement with earlier 
measurements at temperatures above ~ 5 keV. At lower temperatures we see a slightly lower 
normalization than either the 102 or Henry (2000) measurements. 

A particular feature in our cluster sample, that is also evident to some extent in the 
102 and Henry (2000) samples, is a small gap or deficit of clusters in the temperature range 
of ~ 4.5 — 5.5 keV. This is also apparent in Figure 4 and in Figure 6 where we plot cluster 
temperature versus redshift for our sample. In particular, the Suzaku measurements appear 
grouped towards 4-4.5 keV across a range in redshift. Although this may suggest a systematic 
issue in our temperature estimates we have also plotted in Figure 6 the results of our identical 
analysis of Suzaku data on the archived objects A1060, A2218, and A1795 (Table 4), which 
agree very well with pre-existing measurements from other missions, and show no obvious 
biases. This gap is also apparent in the pre-existing temperature measurements that make 
up the rest of our sample. Although the 102 and Henry 2000 cumulative XTF's also show 
a similar feature it is notable that the Markevitch (1998) XTF does not. All samples have 
very similar redshift limits (i.e. z < 0.2 with very few objects beyond z = 0.1). Markevitch 



VmaxiT, L) — — -D^om(^niax) 



(5) 



47r(l + 2;)^£'com(^/)/lim = ^obs(^'^' ^/) 



(6) 
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(1998) derives ASCA temperatures by explicitly removing emission from any cool cores 
(previously ascribed to large cooling flows), whereas 102 and our present work employ a 
similar 2T/1T spectral approach and Henry (2000) used many collated temperatures. It is 
therefore possible that systematic cfi'ccts in temperature measurements are responsible for 
this feature. We do how ever believe that our 1T/2T modeling approach should avoid any 
bias towards underestimating temperatures, since we do detect and reduce the contribution 
of any cool cores. 

It is of course possible that this feature is simply a statistical fluke, due to the flnite 
sample size and L — T scatter. Indeed, there is a clear steepening in the cumulative XTF 
slope (Figure 5) at temperatures immediately below 4.5 keV compared that above 5.5 keV, 
indicating a slight "pile-up" of objects cooler than 4.5 keV. If we naively average over the 
temperature range 2-8 keV then there are approximately 4.5 clusters per 0.5 keV bin. Finding 
zero clusters in the 4.5-5 keV range (see Table 5) is therefore (assuming Poisson statistics) 
signiflcant at about a 2 — o" level. Finding 10 clusters in the 4-4.5 keV range (see Table 5) 
is a positive deviation at about a 2 — 3o" level. Taken together this suggests that while a 
statistical fluctuation cannot be ruled out, it would certainly be well outside a 2-a event. 

It is also possible that this is a manifestation of cosmic variance due to large-scale 
structures within the flnite volume of our sample (which is distinct from a pure statistical 
fluke). This is something that has plagued very early efforts to make cosmological methods 
with galaxy clusters (e.g. Edge et al. 1990). The apparent dearth of objects at 2; ~ 0.05 
(Figures 1 and 6) is likely a consequence of this, although this does not obviously map to 
the deflcit of 4.5-5.5 keV clusters. 



6. Summary and Conclusions 

We have conducted an X-ray galaxy cluster survey of 14 systems using Suzaku. These 
clusters are part of a statistically complete sample of 61 clusters in the local (z < 0.1) Uni- 
verse with fx > 1.0 X 10^^^ erg s~^ cm~^ (0.1-2.4 keV), in the northern celestial hemisphere. 
None of these 14 clusters have been observed previously in pointed X-ray observations, de- 
spite being luminous systems {L^ > 1 x 10^^ erg s~^ (0.1-2.4 keV)). The Suzaku data is 
shown to be of high quality and enables robust imaging and spectral analyses to be made. A 
comparison of spectral analyses of Suzaku data on 3 well known clusters (not in our sample) 
with those from earlier missions indicates no obvious biases either in the Suzaku calibration 
or the analysis approach that we use. 

We present the first results of a spectral analysis of these 14 clusters with the aim of 
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deriving gas temperatures that reflect the true virial temperature. To do this we employed a 
spectral modeling scheme that determines whether a one or two plasma temperature model 
is appropriate. Two temperature models allow us to separate out the contribution of cooling 
cores in 7 of the clusters. Combining the new Suzaku data with preexisting temperature data 
for the flux limited sample we obtain a dataset of 57 objects with spectroscopically conflrmed 
gas temperatures and 4 objects for which pointed X-ray data is either not yet public or has 
not yet been analyzed. For these 4 systems we simply use estimated temperatures based on 
the L — T relation, and note that their inclusion or exclusion has a negligible impact on our 
results. 

We have determined both the L — T relation and the cumulative XTF for the full sample. 
In fltting the L — T relationship with a power law we argue that the often used BCES method 
can produce a power law index that is dependent on the actual measurement errors in L and 
T, owing to the weighting scheme employed. We have therefore also employed a modifled 
flt that takes intrinsic scatter into account in an attempt to mitigate any bias in the 
flt. We flnd that this modifled ^-nd the BCES method are in quite close agreement, with 
power-law indices differing by ~ 4% and normalizations by ~ 0.1 %. The modifled yields 

-t^0.1-2.4keV — J-U 1 

The measured XTF is in close agreement with that of previous studies. Obviously 
there is signiflcant overlap in the actual clusters contained in all such works. However, our 
lower flux limit (by a factor of 2 from, for example, 102) and the fact that 22% of our 
sample has never been studied in the X-ray before, helps reduce this problem. There is 
some evidence that our XTF is slightly lower at temperatures less than ~ 4.5 keV than the 
earlier measurements. This appears to be partially a consequence of a deflcit of clusters in 
the range of 4.5-5.5 keV. We have considered several possible causes for this. These include; 
a systematic bias in our spectral analyses, a statistical fluke, or a consequence of large-scale 
structure and cosmic variance in this very local sample. Since a similar deflcit is also apparent 
in two earlier XTF measurements (102, Henry (2000)) it may be that a statistical or cosmic 
variance effect is the most likely. 

If local galaxy cluster samples are to be used to construct robust calibrations for future 
high-precision cosmological tests then it will be crucial to ensure that wc understand the flnite 
volume effects (e.g. cosmic variance) that may produce systematic biases. One solution will 
be to push measurements to slightly lower mass scales where clusters are more numerous. We 
have demonstrated that Suzaku is an excellent instrument for doing just that, by exploring 
the population of luminous, nearby, clusters that has nonetheless not yet been fully studied. 

We thank Zoltan Haiman and an anonymous referee for helpful comments. This work 
was supported by NASA Suzaku grant NNX06AI40G, and we thank the entire Suzaku team 
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and the staff of the Suzaku Guest Observer Facihty. 
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Fig. 1. — The X-ray luminosity of all BCS clusters is plotted versus redshift. The flux 
and redshift limits of our sample of 61 objects are shown by the dashed and solid lines 
respectively. Notice the apparent deficit of luminous objects at z ~ 0.05. 
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Fig. 2. — Raw images of clusters (background NOT subtracted) from the XISO camera in 
the 0.5-2 keV energy range together with brightness contours. The colour and contours are 
both on a logarithmic scale. The data were binned by 4, and smoothed with a cr = 5 pixel 
Gaussian. 
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Fig. El — Continued 
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Fig. 3. — Spectra of A1800 and Z8276 and IT model fit in the energy range 0.5-7 keV, 
background has aheady been subtracted. According to our F-test, A1800 is a NCCC, while 
Z8276 is CCC Spectra from XISO, XISl, XIS2 and XIS3 are shown in black, red, green and 
blue, respectively. 
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Fig. 4. — Cluster luminosity in the energy range of 0.1-2.4 keV is plotted versus temperature. 
The new Suzaku data is shown together with the data from the rest of our flux-limited sample. 
Both the BCES best fit L — T relation and that of our modified method are shown. 
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Fig. 5. — Cumulative XTF together with previous resuhs. The thick black solid line is our 
present result, the red dotted line is from 102, the blue short-dashed line is from Markevitch 
(1998), the cyan dash-dotted line is from Henry (2000). The shaded (yellow) region indicates 
the 1-0" poisson error on our XTF measurement. 
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Fig. 6. — Cluster temperature is plotted versus redshift for our sample. Suzaku cluster 
measurements are shown in red. The three objects used to compare our X-ray temperature 
analysis to that of earlier missions (Table 4) are plotted as filled circles. 
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Table 1: Suzaku observations of 14 galaxy clusters 



Cluster Name 


Sequence Number 


(RA,DEC)(deg) 


Exposure Time (ks) 


Processing Version 


UGC03957 


801072010 


(115.2386, 55.4389) 


12.4 


1.2.2.3 


A2572b 


801073010 


(349.6122, 18.7237) 


22.4 


1.2.2.3 


Z8852 


801074010 


(347.6298, 7.5832) 


8.6 


1.2.2.5 


II Zw 108 


801075010 


(318.5265, 2.5682) 


11.4 


1.2.2.3 


A2665 


801076010 


(357.7129, 6.1379) 


10.5 


1.2.2.5 


A1831 


801077010 


(209.8017, 27.9785) 


8.9 


1.2.2.3 


A1800 


801078010 


(207.3647, 28.1066) 


18.5 


1.2.2.3 


Z8276 


801079010 


(266.0553, 32.9813) 


14.9 


1.2.2.3 


A2495 


801080020 


(342.5710, 10.9166) 


11.9 


1.4.2.9 


A2061 


801081010 


(230.3133, 30.6555) 


15.3 


1.2.2.3 


A2249 


801082010 


(257.4518, 34.4734) 


22.8 


1.3.2.6 


Z235 


801083010 


(10.9674, 24.4018) 


12.9 


1.2.2.3 


A272 


801084010 


(28.7653, 33.8978) 


14.1 


1.2.2.3 


A566 


801085010 


(106.0943, 63.2750) 


8.79 


1.3.2.8 
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Table 2: Spectral Fitting Result (Metallicity fixed) 



Cluster Name Redshift nn IT Model 2T Model 









^-| n20 
(10 


cm j 


Ki (keV j 


X 


2 / J „ f 

/d.o.f 


kih[kc\) 


fcTKkeV) 


X 


2 / J „ J" 
/d.o.f 




0- 


no A 1 

.0341 


/I 

4, 


.35 




1495- 


.19/887 


Q ^7 + 0.11 


1, 


Q^+0.07 
•«4_o.03 


1457 


.29/885 


AZO (ZD 


n 

U, 


.Uooy 


A 

4, 


.iU 


70+0.04 
^•'"^-0.04 


iODD, 


.49/00/ 




1, 


0^+0.08 


i4z4 


.yo/ooo 


ZjooOZ 


U, 


.U4UU 


A 

4, 


. / o 


''■51^-0.11 


1 ni Q 
iUio, 


1 1 IRRK 


q 70+0.78 
^•'^-0.96 


3 


17+1.64 
■-^'-0.82 


1 ni Q 
iUio 


.iZ/DDo 


II Zw 108 


0, 


.0483 


5, 


,85 


3.90l°:i° 


1043 


.86/887 


o nv+O.ll 
-0.11 


1, 


1 q+4.08 

•-•^"^-0.61 


1041, 


.17/885 


A2665 


0, 


.0562 


5, 


.82 


3.98l°ii 


1038, 


.76/665 


14 Q1+9.97 


2, 


QO+0.21 


1015, 


.34/663 


A1831 


0, 


.0612 


1, 


.37 


•J-'J'J-o.oe 


1286, 


.31/887 


4 22+0^^ 

^•■^■^-0.14 


1, 


14+0.15 

•1^-0.06 


1178, 


.48/885 


A1800 


0, 


.0748 


1, 


.17 




1196, 


.85/887 


4 iq+0.33 
^•1^-0.46 


3 


07+2.86 
•O'-1.50 


1196 


.85/885 


Z8276 


0, 


.0757 


3, 


.88 


O ro+0.05 

^•'J°-0.05 


1252, 


.13/887 


4 09+° °^ 


1, 


4t^+0.08 
•^^-0.03 


1176 


.05/885 


A2495 


0, 


.0768 


4, 


.82 


^■-'-'^-0.09 


1243, 


.78/665 


o 71+4.33 
O. / J-_2 27 


3 


71+0.55 
• ' -L-0.13 


1241 


.32/663 


A2061 


0, 


.0777 


1, 


.67 


4 sn+°-i° 

^•"^'^-O.lO 


1035 


.92/887 


4 40+0.12 

0.12 


0, 


71+0.14 
• ' -'--0.14 


1024, 


.50/885 


A2249 


0, 


.0802 


2, 


.22 


^ 4fi+0.16 
0.16 


879, 


.93/665 


12.10l|;95 


4, 


f;9+0.42 
■"^-0.22 


875, 


.49/663 


Z235 


0, 


.0830 


3, 


.73 


q c-9+0.06 


1451, 


.43/887 


q QQ+O-io 

«>»y_o.o6 


1 


1 c;+0.07 
•-'-"-'-0.02 


1343, 


.01/885 


A272 


0, 


.0872 


4, 


.99 


o 26+0-07 


1125, 


.81/887 


99+0.16 


1, 


71+0.10 
• ' -'--0.06 


1104 


.59/885 


A566 


0, 


.0980 


5, 


.03 


o 7Q+0.13 


1098, 


.36/665 


43.47 


3 


■^O-O.IO 


1096 


.52/663 



Notes — Errors in the tables are all 90% error. 
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Table 3: Spectral Fitting Result (Metallicity as free parameter) 



Cluster Name 


IT Model 


2T Model 




A;T(keV) 


Z{Zq) xVd.o.f 


A;T,,(keV) kTi(keY) Z{Zq) 


xVd.o.f 



UGC03957 

A2572b 

Z8852 

II Zw 108 

A2665 

A1831 

A1800 

Z8276 

A2495 

A2061 

A2249 

Z235 

A272 

A566 



-^•O" 0.04 



0.59^ 



0.05 
-0.05 
0.11 
0.11 
0.10 
0.10 
+0.11 
0.11 
+0.06 
0.06 
+0.09 
-0.09 



2.72 
3.50 
3.93 
3.96 
3.64 
4.13 

+0.09 
0.09 
+0.11 
-0.10 
+0.16 
-0.16 
+0.05 
-0.05 
+0.07 
-0.07 
+0.13 
-0.12 



4.10 

4.31 
5.52 
3.51 
3.26 
3.79 



hO.03 
-0.03 

"-'•^0-0.02 

n. qo + 0.06 

+0.04 
0.04 
0.04 
0.04 

+0.03 
-0.03 
0.03 
0.03 
+0.03 
0.03 
0.04 
0.03 



0.45 
0.35 
0.31 
0.36 
0.50 
0.35 

"-'•^^-0.03 
+0.04 
-0.04 
+0.03 
-0.03 
+0.03 
0.03 
+0.05 
0.04 



0.20 
0.45 
0.31 
0.33 



1228.23/886 
1560.63/886 
1012.76/664 
1009.83/886 
1035.00/664 
1270.34/886 
1189.41/886 
1109.93/886 
1237.65/664 
1027.97/886 
862.40/664 
1368.40/886 
1125.61/886 
1096.84/664 



(-4+0.31 


z, 


,55_o 


.19 
.35 


0-59_o,03 


•*''^-0.26 


1 

1, 


qc+O, 


.07 
.12 


n 90+0.02 

"•'^''-0.03 


q «q+0.82 


3, 




.65 
.82 


n 00+0,06 
'-'•'^^-0,05 


4.03^^:}^ 


0, 


.68^°: 


.22 
.16 




18.30^ 


3, 


.071°: 


.17 
.15 




4 j^g+0.20 


1, 




.15 
.07 


Q 01 +0.03 
^••'^-0.03 


4 iq+0.82 
^•■'-^-0.82 


3, 




.29 
.82 


S6+° °3 


4.08lH^ 


1, 


.711°: 


,09 
.10 


U.410 Q Q3 


11.98ltf? 


3, 


.601°: 


.13 
.09 


n qo+0,04 
'^•'^^-0,04 


4.44l°ii 


0, 


,711°: 


.19 
.15 


"•■^'^-0.03 




5, 


,431°: 


.82 
.82 


'-'•"^'-'-0.04 


3.801HI 


1, 


,141°: 


.05 
.05 


"•^"-0.04 


4 Qo+0.16 


1, 


,671°: 


.11 
.13 


26+° °^ 
"•■^"-0.03 


31.16 


3, 


,671°: 


.13 
,10 


"-'•^^-0.04 



1228.24/884 

1409.27/884 

1012.77/662 

1003.58/884 

1000.71/662 

1161.77/884 

1189.43/884 

1081.53/884 

1229.03/662 

1018.67/884 

862.42/662 

1292.94/884 

1101.16/884 

1094.36/662 
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Table 4: Calibration clusters 



Cluster Name 


Tnew(keV) 


^old(keV) 


Reference 


A1060 


3.30 


3.1-3.9 


3,4,6,7,8,13,14,15 


A2218 


6.90 


6.7-7.2 


1,3,5,9,11,12 


A1795 


5.73 


5.3-7.8 


2,3,8,10,16,17 



Reference— (l)Allen 1998; (2)Allen et al. 2001;(3)David et al. 1993; (4)Furusho et al. 2001; (5)Govoni et 

al. 2004; (6)Hayakawa ct al. 2004; (7)Hayakawa et al. 2006; (8)102; (9)Macliacek et al. 2002; (lO)Markevitch 
1998; (ll)McHardy et al. 1990; (12)Mushotzky & Loewenstein 1997; (13)Sato et al. 2007;(14)Tamura et al. 
1996; (15)Tamura et al. 2000; (16)Taniura et al. 2001; (17)Viklilimn et al. 2005 
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Table 5: Complete Local Flux Limited Cluster Sample 



Cluster Name 


Redshift 


(keV) 


Flux 

(lO-l^erg cm-2 s"^) 


Luminosity 
(10''4/i-2erg s-l) 


(/i.-3Mpc-3) 


Instrument 


Reference 


Virgo 


0.0036 


-0.52 


1821.1 ±38.8 


0.254 ± 0.005 


1.23E+07 


ASCA 


11 


RXJ0419.6+0225 


0.0123 


1 4+0.06 
^•^-0.06 


45.7 ±5.4 


0.075 ± 0.009 


4.03E+05 


ASCA 


8 


A262 


0.0163 


9 9rr+0.06 


48.6 ± 4.8 


0.141 ±0.014 


2.20E+06 


ASCA 


7 


RXJ0123.6+3315 


0.0164 


1 97+0.04 
r-^'_0.04 


17.1 ±4.0 


0.050 ±0.012 


2.82E+05 


ASCA 


1 


Z4803 


0.0200 


1 oc+0.02 
^•°°-0.02 


17.7 ±2.2 


0.078 ± 0.010 


l.lOE+06 


ASCA 


5 


A1367 


0.0214 


Q cc+O.OS 
"^•°°-0.08 


81.4 ±4.9 


0.408 ± 0.024 


l.lOE+07 


ASCA 


7 


RXJ0751.3+5012 


0.0220 


1.98= 


10.1 ± 1.8 


0.054 ±0.010 


1.40E+06 






A1656 


0.0231 


8.071°;^? 


319.2 ±8.3 


1.863 ± 0.048 


3.50E+07 


ASCA 


7 


A400 


0.0238 


2.43i°;l| 


16.2 ± 3.3 


0.100 ±0.021 


2.90E+06 


ASCA 


7 


MKW8 


0.0276 


n QQ+0.23 

'-'•^^-0.22 


18.4 ±2.7 


0.154 ±0.023 


8.50E+06 


ASCA 


7 


RXJ1715.3+5725 


0.0280 


2 2fi+''--'-'' 
^•^"-0.06 


12.7 ±0.9 


0.110 ±0.008 


2.24E+06 


ASCA 


6 


RXJ0341.3+1524 


0.0290 


2.96= 


15.5 ±3.1 


0.142 ±0.028 


5.83E+06 






A2199 


0.0299 


4 28+"-^*' 
*"'°-0.10 


96.8 ± 3.5 


0.955 ± 0.034 


1.97E+07 


ASCA 


7 


A2634 


0.0309 


3 4ci+016 
•3-^o_0.i6 


23.1 ±2.5 


0.243 ± 0.026 


l.OOE+07 


ASCA 


7 


AWM4 


0.0318 


9 c:fl+0 04 


10.2 ± 1.5 


0.114 ±0.017 


3.49E+06 


XMM 


10 


RXJ1733.0+4345 


0.0330 


1 tro+O.OS 


11.6 ± 1.1 


0.140 ±0.013 


5.42E+05 


Chandra 


4 


UGC03957 


0.0341 


o .7+0.07 
■^•^'-0.12 


18.6 ±2.6 


0.239 ± 0.033 


1.02E+07 


Suzaku 




2A0335+096 


0.0349 


„ ^ j+0 09 
^■(^-0.08 


80.5 ± 6.5 


1.073 ± 0.087 


1.20E+07 


ASCA 


7 


A2147 


0.0353 


4 34+0.12 


53.2 ±3.7 


0.732 ± 0.051 


2.04E+07 


ASCA 


7 


A2052 


0.0353 


Q 19+0.10 


47.1 ±3.7 


0.650 ± 0.051 


7.06E+06 


ASCA 


7 


A2063 


0.0355 


3.56l°;l^ 


37.5 ±3.4 


0.523 ±0.048 


1.12E+07 


ASCA 


7 


A2151a 


0.0370 


3.80t°;I° 


15.0 ± 1.9 


0.227 ± 0.028 


1.38E+07 


Einstein 


2 


A576 


0.0381 


Q QQ+0.16 

3-»3_() 15 


22.4 ±3.5 


0.360 ± 0.056 


1.42E+07 


ASCA 


7 


A2572b 


0.0389 


Q oa+0.08 
3-^°-0.12 


10.0 ± 1.9 


0.168 ±0.031 


8.40E+06 


Suzaku 




A76 


0.0395 




11.0 ±2.4 


0.192 ± 0.041 


5.17E+05 


Einstein 


2 


A2657 


0.0400 


n c:q+0.12 
O.JO_Q 


23.4 ± 2.9 


0.415 ±0.052 


1.08E+07 


ASCA 


7 


Z8852 


0.0400 


cn+0.07 


13.8 ±2.1 


0.245 ± 0.037 


1.05E+07 


Suzaku 




A2107 


0.0411 


QO+0.10 

•J-aj-o.io 


15.2 ±2.4 


0.285 ± 0.045 


1.54E+07 


ASCA 


9 


A2589 


0.0416 


00+0.13 

•J-JO_0.13 


25.3 ± 2.8 


0.487 ± 0.055 


9.34E+06 


ASCA 


7 


HCG94 


0.0422 


on+0.17 
O.OU_(, ig 


17.2 ± 0.03 


0.341 ± 0.006 


8.60E+06 


ASCA 


7 


A2593 


0.0428 


•^•^"-0.90 


14.8 ±2.1 


0.302 ±0.044 


6.90E+06 


Einstein 


2 


A168 


0.0448 


r, co + 0.09 


11.0 ± 1.8 


0.247 ± 0.041 


3.35E+06 


Chandra 


12 


MKW3S 


0.0453 


4r+0.13 


30.3 ± 3.3 


0.694 ± 0.075 


l.OOE+07 


ASCA 


7 


IIZwl08 


0.0483 


, 0(1+0.06 


21.9 ±2.8 


0.569 ± 0.073 


1.50E+07 


Suzaku 




A376 


0.0488 


, 7C+0.05 


13.6 ±2.0 


0.360 ± 0.053 


1.32E+07 


ASCA 


3 


A193 


0.0491 




15.4 ±3.0 


0.414 ±0.080 


1.87E+07 


EXOSAT 


2 


A2665 


0.0562 


3.98iro? 


14.1 ±2.2 


0.499 ± 0.076 


1.60E+07 


Suzaku 
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Cluster Name 


Redshift 


(keV) 


Flux 

(10~^^erg cm~^ s~^) 


Luminosity 
(lO'l'i/j-aerg s"!) 


(/j-3Mpc-3) 


Instrument 


Reference 


A2626 


0.0565 


2-90+^:g'(j 


14.3 ± 2.0 


0.515 ±0.071 


5.45E+06 


Einstein 


2 


A2256 


0.0581 


6-83+0 21 


49.5 ± 2.2 


1.867 ± 0.082 


3.46E+07 


ASCA 


7 


A1831 


0.0612 




11.9 ±1.6 


0.505 ± 0.067 


1.89E+07 


Suzaku 




A 1795 


0.0622 


6.17i«-^6 


68.1 ±4.0 


2.976 ± 0.173 


3.37E+07 


ASCA 


7 


A1767 


0.0701 


4.10tl:?° 


11.7 ± 1.4 


0.656 ± 0.076 


1.74E+07 


Einstein 


2 


A1589 


0.0718 


5.44^ 


10.8 ±1.6 


0.634 ± 0.094 


3.09E+07 






A399 


0.0722 


6.46l°i 


28.8 ± 3.5 


1.709 ± 0.208 


3.42E+07 


ASCA 


7 


A2065 


0.0723 


r 07+0.34 

"-'•"J ' _o 30 


22.2 ±2.8 


1.317 ±0.166 


3.04E+07 


ASCA 


7 


A1775 


0.0724 




13.0 ± 1.9 


0.782 ± 0.112 


1.22E+07 


ASCA 


7 


A401 


0.0739 


7 iqO.28 


42.8 ±3.9 


2.607 ±0.237 


3.48E+07 


ASCA 


7 


A1800 


0.0748 


4 14O.O6 
—0.06 


12.8 ±1.8 


0.822 ±0.115 


1.79E+07 


SuzaJju 




Z5029 


0.0750 


7.52*^' 


22.0 ± 2.4 


1.404 ±0.155 


3.49E+07 






Z8276 


0.0757 




16.4 ±1.7 


1.604 ±0.111 


1.73E+07 


Suzaku 




A2029 


0.0766 


7 03+0:04 


61.3 ±4.3 


4.068 ± 0.288 


3.50E+07 


ASCA 


7 


A2495 


0.0768 


4 120-05 


11.8±2.1 


0.792 ± 0.141 


1.77E+07 


Suzaku 




A2061 


0.0777 




15.3 ±2.2 


1.059 ±0.153 


2.14E+07 


Suzaku 




A2249 


0.0802 




14.4 ± 1.6 


1.059 ±0.115 


3.10E+07 


Suzaku 




A2255 


0.0809 


■^-"^-0.26 


17.8 ±8.9 


1.330 ± 0.665 


3.30E+07 


ASCA 


7 


Z235 


0.0830 


4 31 +0.06 


10.9 ±1.8 


0.856 ±0.140 


2.00E+07 


Suzaku 




A272 


0.0872 


gq+0.10 


10.1 ± 2.1 


0.881 ± 0.186 


1.61E+07 


Suzaku 




A478 


0.0882 


6 qi+"-'"' 


39.9 ±5.1 


3.465 ± 0.441 


3.47E+07 


ASCA 


7 


A2142 


0.0894 


4f; + 0.53 

0-46_o 4y 


61.5 ±3.9 


5.589 ± 0.357 


3.50E+07 


ASCA 


7 


A2244 


0.0970 


t- -7+0.61 
=5-' '-0.44 


23.4 ± 2.0 


2.542 ±0.215 


3.25E+07 


ASCA 


7 


A566 


0.0980 


7q+0.13 
"^-'^-0.12 


11.6 ±2.2 


1.284 ± 0.239 


1.35E+07 


SuzaJfu 





Notes — All the errors are 1 a error. The 90% temperature error in literature is converted to 1 cr error 
assuming it's Gaussian error. 
"Temperature with superscript e is temperature estimated from our L — T relation. 
REFERENCE— (l)Buote & Fabian 1998 ; (2)David et al. 1993; (3)Dupke & Bregman 2005; (4)Fukazawa et 
al. 2006;(5) Horner 2001; (6)Hwang ct al. 1999; (7)102; (8)Kawaharada et al. 2003; (9)Novicki et al. 2002, 
(lO)O'Sullivan et al. 2005; (ll)Shibata et al. 2001; (12)Yang et al. 2004. 
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Table 6: L — T Relation Fitting Results 



Fitting Method 


Lq.i 


-2.4 — T 




L 


Bol — T 






a 


P 


'^intr 


a 


p 


^intr 




42.31 


2.45 


0.18 


42.31 


2.95 


0.18 


BCES 


42.26 ±0.12 


2.55 ±0.19 


0.18 


42.20 ±0.14 


3.16 ±0.22 


0.18 



